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Abstract- In  the  transmission  CT  imaging  system  with  scintilla- 
tor-CCD  detector  using  synchrotron  monochromatic  x-ray  as  sheet 
beam,  improvement  of  reconstructed  image  quality  is  investi¬ 
gated.  The  properties  of  experimental  facilities  cause  the  inci¬ 
dent  beam  intensity  to  be  spatially  inhomogeneous,  and  to  tem¬ 
porally  attenuate.  In  order  to  achieve  more  reliable  reconstruc¬ 
tion,  we  propose  the  data-correction  and  estimation  method  based 
on  the  properties  of  the  facilities  and  the  principle  of  CT  recon¬ 
struction,  and  confirm  its  effectiveness  using  a  physical  phan¬ 
tom. 

I.  Introduction 

Synchrotron  radiation  (SR)  attracts  researchers  in  biomedical 
imaging  field  because  of  its  excellent  properties  such  as  broad 
and  continuous  spectrum,  high  brilliance,  and  high  directional¬ 
ity  From  the  properties,  SR  provides  much  higher-quality  mono¬ 
chromatic  x-ray  beam  with  high  brilliance  and  directionality  than 
the  conventional  x-ray  tube,  and  then  releases  the  difficulty  in 
imaging  by  x-ray  tube,  e.g .,  beam  hardening  effects  causing  deg¬ 
radation  of  image  qualities  [1-14].  To  obtain  high  contrast  CT 
images,  there  is  introduced  a  highly  sensitive  imaging  system 
using  scintillator-CCD  detector  system  by  incident  sheet  beam, 
which  drastically  reduces  data-acquisition  time  compared  with 
pencil  beam,  while  sheet-beam-used  imaging  slightly  degrades 
image  qualities  due  to  scattered  x-rays  [7,  8,  12,  13,  14].  In  spite 
of  using  both  the  high-quality  x-ray  source  and  the  highly  sen¬ 
sitive  imaging  system,  the  reconstructed  image  quality  is  de¬ 
graded  because  of  the  following  reasons:  First,  the  intensity 
profile  of  incident  sheet  beam  represents  spatial  inhomogeneity. 
Second,  the  intensity  of  incident  beam  temporally  attenuates  in 
the  course  of  measurements  due  to  the  properties  of  experimen¬ 
tal  facilities.  The  fact  reveals  that  high  quality  SR- imaging  re¬ 
quires  not  only  highly  sensitive  imaging  system  but  also  reli¬ 
able  data  estimation  based  on  the  properties  of  the  experimental 
facilities.  In  this  research,  we  propose  the  data-estimation  and 
correction  method  for  SR-imaging  with  incident  sheet  beam,  and 
demonstrate  its  effectiveness  using  a  physical  phantom. 

II.  Transmission  x-ray  ct  imaging  system 

BY  INCIDENT  SHEET  BEAM 

Fig.  1  shows  the  constitution  of  transmission  CT  imaging 
system  by  incident  x-ray  sheet  beam  [7, 14].  Incident  white  x-ray 
beam  from  the  accelerator  is  monochromated  with  the  Si  (1 1 1) 


double-crystal  Bragg-Bragg  monochromator.  Then,  the  object 
on  the  rotational  scanning  stage  is  irradiated  with  the  beam  col¬ 
limated  to  2x40  mm2  rectangular  intensity  profile  using  the  tan¬ 
talum  slit.  The  transmitted  x-ray  is  converted  to  visible  light  by 
the  scintillating  screen  for  medical  purposes  (HR-6,  FUJI  FILM), 
which  is  detected  with  the  CCD  camera  (CCD  15-11,  Astro  Cam) 
after  being  reflected  by  the  mirror  and  focused  by  the  lens.  The 
CCD  and  the  rotational  stage  are  controlled  by  PC.  The  optical 
system  from  the  scintillating  screen  to  the  CCD  is  covered  with 
black  sheet  to  avoid  strayed  light,  and  the  CCD  is  set  at  right 
angle  against  the  incident  beam  to  avoid  the  direct  x-ray  irradia¬ 
tion,  and  detects  signal  light  as  2-D  data  with  1 6-bit  gray  scale, 
whose  pixel  size  and  image  size  are  181x181  pm2,  and  1024x256 
pixels,  respectively.  The  system  was  constructed  at  bending 
magnet  beam  line  NE-5A  of  the  Tristan  Accumulation  Ring  (6.5 
GeV,  10-40  mA)  in  KEK,  Japan.  Fig.  2  shows  an  example  of  mea¬ 
sured  2-D  data  without  an  object  and  its  vertical  and  horizontal 
profiles.  Projection  data  was  generated  from  the  2-D  image  data 
from  the  CCD  as  follows:  First,  irradiated  region  is  extracted  by 
the  thresholding  method,  and  next  the  datum  in  a  bin  is  gener¬ 
ated  by  summing  pixel-values  of  adjacent  12  and  10  pixels  in 
vertical  and  horizontal  directions,  respectively,  for  reducing  sta¬ 
tistical  noise.  Thereby,  projection  data  were  obtained. 

Anti-reflection 


Fig.  1.  Schematic  of  tansmission  x-ray  computed  tomography  imag¬ 
ing  system  using  synchrotron  monochromatic  x-ray  sheet  beam. 

III.  Correction  for  spatial 

INHOMOGENEITY  IN  INCIDENT  BEAM  INTENSITY 

Fig.  3  shows  a  projection  without  an  object  obtained  from  the 
data  of  Fig.  2.  One  can  see  the  spatial  inhomogeneity  of  incident 
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Fig.  2.  Example  of  measured  incident  beam  and  its  vertical  and  hori¬ 
zontal  profiles 


Position  [pixel] 


Fig.  3.  Example  of  spatially  inhomogeneous  indident  intensity 
obtained  from  the  measured  data 

beam.  The  causes  of  the  inhomogeneity  are  as  follows:  (1)  inho¬ 
mogeneity  of  sensitivity  among  pixels  of  CCD  camera,  (2)  spatial 
inhomogeneity  of  incident  white  beam  intensity  itself  from  an 
accelerator,  and  (3)  monochromator  deformed  by  heat  due  to 
incident  white  beam.  For  (1),  one  can  easily  correct  the  effect  by 
employing  projections  obtained  from  a  water-filled  cylindrical 
physical  phantom.  We  do  not  consider  (1)  because  the  effective 
correction  method  has  been  established.  For  (2),  it  is  very  diffi¬ 
cult  to  obtain  homogeneous  incident  beam  profile  in  spite  of 
carefully  adjusting  it  before  setting  up  the  system.  For  (3),  al¬ 
though  the  monochromator  is  water-cooled  to  reduce  the  heat 
effect,  it  is  impossible  to  suppress  the  deformation  because  of 
continuous  irradiation  of  white  beam  with  intense  intensity.  (2) 
and  (3)  are  proper  to  synchrotron  radiation  facilities. 

Since  the  homogeneities  in  incident  beam  intensity  caused 
by  (2)  and  (3)  depend  on  the  status  of  experimental  facilities 
when  imaging,  correction  method  must  reflect  the  status.  So,  we 
paid  attention  to  an  intensity  profile  obtained  just  before  imag¬ 
ing  an  object.  Let  z(r),  t(r,  6 ),  and  t  '(r,  6)  be  an  intensity  profile  of 
incident  beam  obtained  just  before  measurements,  a  measured 
transmission  profile  at  projection  angle  6 ,  and  its  corrected  trans¬ 
mission  profile,  respectively,  where  r  and  0are  defined  as  shown 
in  Fig.  4.  Assuming  that  the  spatial  inhomogeneity  does  not 
change  in  the  course  of  measurements,  one  can  estimate  the 
transmission  profile  t  '(r,  0)  obtained  when  imaging  the  object 
with  spatially  homogeneous  sheet-beam  by  normalizing  the  mea¬ 


sured  transmission  profile  t(r ,  0)  with  an  intensity  profile  i(r)  as 
t  \r ,  0)  =  t(r ,  0)  iQ  /  i(r ),  (1) 

where  zQ  is  a  spatial  average  value  of  z(r). 

IV.  Correction  for  temporal  attenuation 

IN  INCIDENT  BEAM  INTENSITY 

With  a  decrease  in  the  number  of  electrons  in  the  accelerator 
ring  the  synchrotron  radiation  beam  intensity  gradually  decreases 
proportionally  to  the  electron  number.  Fig.  5  shows  the  temporal 
attenuation  of  projection  data  obtained  from  a  physical  phan¬ 
tom.  One  way  in  which  one  detects  the  intensity  attenuation  is  a 
ring  current  value  proportional  to  the  number  of  electrons  in  the 
ring.  In  fact,  so  far  we  have  employed  the  ring  current  data  to 
correct  the  intensity  attenuation.  However,  although  the  ring 
current  is  highly  correlated  with  the  beam  intensity,  there  does 
not  exist  the  strictly  quantitative  relationship  between  them. 
Thus,  it  is  insufficient  for  a  correction  parameter.  Another  way  is 
to  monitor  the  attenuation  with  an  ion  chamber  set  on  the  upper- 
stream  side  than  the  object.  However,  although  the  method  is 
effective  for  measurements  with  pencil  beam,  it  is  not  suitable  to 
measurements  with  sheet  beam  because  an  ion  chamber  can  not 
acquire  a  spatial  information  of  the  beam. 

First,  we  pay  attention  to  the  fact  that  total  X-ray  quantity 
absorbed  in  the  object  is  constant  when  imaging  with  the  spa¬ 
tially  homogeneous  incident  sheet  beam  independently  of  pro¬ 
jection  angle  6.  In  fact,  letting  /(x,  y)  be  the  spatial  distribution  of 
the  linear  attenuation  coefficients,  the  projection  for  projection 
angle  6  is  given  as 
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Fig.  5.  Example  of  temporal  attenuation  of  transmitted  x-ray 
intensity 


g(r,6) 


»+oo 

:  J  f(r  cosO  -  s  sin6,r  sinO  -  s  cos6)ds 


(2) 


Accordingly,  the  above-mentioned  fact  can  be  easily  proved  as 
follows: 


g(r9e) 


-  J  J  f(rcos0-ssin0,rsin0-scos0)dsdr 

f{x>y)dxdy  =  const. ,  (3) 


where  the  integral  interval  was  set  to  (-©o5  oo)  because  f(x,  y)  is  a 
spatially  localized  function,  and  the  relationship  between  coor¬ 
dinate  systems,  (x,  y)  and  (x,  t ),  was  used. 

The  relationship  among  the  transmission  profile  corrected 
for  the  spatial  inhomogeneity  t'(r ,  0),  the  average  intensity  iQ, 
and  the  projection  g(r,  6)  is  given  as 

log t'(r,e)  =  logi0-g(r,0) .  (4) 

Here,  note  that  the  corrected  intensity  t'(r ,  6)  is  assumed  to  be 
obtained  using  spatially  homogeneous  incident  sheet  beam.  In¬ 
tegrating  with  respect  to  r,  and  then  arranging, 

L(0)  (log  r u ■ log '°  ^r=~\Rm  9^r  ’  © 

where  integral  interval  R(6 ),  which  depends  on  projection  angle 
6 ,  is  introduced  as  shown  in  Fig.  4  to  reduce  to  the  utmost  the 
influence  of  the  inhomogeneity  of  the  region  outside  the  object 
in  the  intensity  upon  the  correction.  If  incident  intensity  does 
not  temporally  decrease,  the  left  hand  side  of  Eq.  (5)  is  constant 
independently  of  projection  angle  0from  the  above  fact.  There¬ 
fore,  one  should  correct  the  temporal  attenuation  such  that  j  log 
0)/iQ)  dr  for  each  6  be  constant  by  normalizing  it  to  the 
value  at  projection  angle  6=  0, 1  log  0 )/z0)  dr.  Hence,  the 

transmission  intensity  is  corrected  as  follows: 


\ogt"(r,0)  =  log(t'(r,0)  /  i0) 

XL(0)  logf'M)  /  *o  Yr  1  Jjj(e)  lQg (*'(r ,  /  h  )dr  +  log  i0 .  (6) 

Fig.  6  shows  a  result  in  which  the  transmission  intensities  of  Fig. 
5  are  corrected  with  the  proposed  method.  One  can  see  that  the 
spatial  and  temporal  inhomogeneities  are  satisfactorily  corrected. 
Finally,  the  projection  data  is  estimated  as  follows: 

g'(r,e)  =  \ogi0-\ogt"(r,e).  (7) 
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Fig.  6.  Profiles  of  transmission  x-ray  intensity  after 
correction 


V.  Experimental  result 

Fig.  7  shows  the  35-mmc|)  cylindrical  physical  phantom  made 
of  acrylic  which  has  three  channels  with  diameter  of  3, 2, 1 .5,  and 
1  mm,  respectively,  i. e. ,  totally  12  channels.  The  incident  white  x- 
ray  is  monochromated  to  34.6  keV.  The  exposure  time  of  CCD 
detector  is  20  ms  so  as  to  effectively  use  its  dynamic  range.  The 
phantom  is  imaged  over  180  degrees  with  rotational  step  of  1 
degree,  and  the  image  is  reconstructed  using  the  filtered 
backprojection  method  with  the  Shepp-Logan  filter.  Figs.  7  (b), 
and  (c)  show  the  reconstructed  images  with  the  ring  current, 
and  the  proposed  corrections,  respectively.  From  the  images, 
one  can  see  that  the  image  with  the  proposed  correction  is  re¬ 
constructed  more  smoothly  than  that  with  the  conventional 
correction  based  on  the  ring  current  value.  To  confirm  the  fact, 
the  profile  of  the  reconstructed  images  on  the  line  passing 
through  the  centers  of  three  channels  of  1  mmc|)  is  shown  in  Fig. 
8,  where  the  large  peak  shown  in  the  right  region  of  the  profile 
for  the  ring  current  correction  is  truncated  at  level  255  for  de¬ 
tailed  comparison  of  image  contrast.  From  the  figure,  one  can 
again  see  that  the  acrylic  region  of  the  image  with  the  proposed 
correction  is  more  smooth  than  that  with  the  ring  current  correc¬ 
tion.  Furthermore,  one  can  recognize  that  the  edges  of  the  image 
with  the  proposed  correction  is  sharper  than  that  with  the  ring 
current  correction. 


Proceedings  -  23rd  Annual  Conference  -  IEEE/EMBS  Oct.25-28,  2001,  Istanbul,  TURKEY 


3-3mm<|) 
3-2mm{() 
3-1.5mm<|> 
3-1mm<|> 

35mrru!> 


Fig.  7.  Comparison  of  reconstructed  images,  (a)  Physical  phan¬ 
tom,  (b)  Reconstructed  image  of  the  physical  phantom  with 
correction  by  ring  current,  (c)  Reconstructed  image  of  the  physi¬ 
cal  phantom  with  correction  of  spatial  inhomogeniety  and  tem¬ 
poral  attenuation. 
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Fig.  8.  Horizontal  profiles  of  the  reconstructed  images 
of  the  physical  phantom 


V.  Concluding  remarks 

In  synchrotron  monochromatic  x-ray  sheet  beam  used  trans¬ 
mission  CT  imaging  system  with  the  scintillator-CCD  detector, 
we  proposed  and  investigated  the  corrections  for  spatial  homo¬ 
geneity  in  incident  beam  intensity  by  using  the  incident  beam 
profile  obtained  just  before  imaging  objects,  and  for  temporal 
attenuation  in  incident  beam  intensity  by  using  the  constant 
value  based  on  the  CT  principle.  As  a  result,  it  is  indispensable 
to  correct  and  then  estimate  the  projection  data  based  on  the 
properties  of  the  experimental  facilities  and  the  CT  reconstruc¬ 
tion  principle  in  order  to  reliably  reconstruct  images  with  the 
excellent  abilities  of  synchrotron  radiation  as  the  light  source. 

However,  Fig.  6  shows  that  the  regions  outside  the  object  in 


the  1 80th  estimated  projection,  where  x-rays  directly  impinge  on 
the  detector  without  passing  through  the  object,  are  different  in 
shape  from  those  at  other  projections.  This  fact  represents  that 
the  assumption  that  the  incident  intensity  profile  does  not 
change  throughout  the  measurement  does  not  hold.  Although 
we  introduced  the  integral  interval  R(6)  to  reduce  this  difficulty, 
it  may  be  insufficient  for  more  reliable  reconstruction.  It  will  be 
effective  to  correct  projections  using  two  profiles,  i.e.,  the  inten¬ 
sity  profiles  without  an  object  obtained  just  before  and  after  the 
measurement,  while  we  used  only  one  profile  obtained  just  be¬ 
fore  the  measurement. 
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